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Abstract 

In this paper, we present a global state feedback tracking controller for underactuated surface marine 
vessels. This controller is based on saturated control inputs and, under an assumption on the reference 
trajectory, the closed-loop system is globally asymptotically stable (GAS). It has been designed using 
a 3 Degree of Freedom benchmark vessel model used in marine engineering. The main feature of our 
controller is the boundedness of the control inputs, which is an essential consideration in real life. In 
absence of velocity measurements, the controller works and remains stable with observers and can be 
used as an output feedback controller. Simulation results demonstrate the effectiveness of this method. 

Index Terms 

Global tracking, bounded feedback, Lyapunov function, underactuated surface marine vessels. 

I. Introduction 

Precise tracking control of surface marine vessels (ships and boats) is often required in critical 
operations such as support around off-shore oil rigs UJ- This problem is of particular interest as 
marine vessels are often underactuated, i.e. the number of independent actuators is less than the 
degrees of freedom (DOF) to be controlled. In this paper, we consider the problem of tracking 
control of a 3-DOF vessel model (surge, sway and yaw 0), working under two independent 
actuators capable of generating surge force and yaw moment only. It has been shown in (3), 
0, AS that under Brockett's necessary condition 0, stabilization of this system is impossible 
with continuous or discontinuous time-invariant state feedback. This can be seen in j6l where the 
authors developed a continuous time-invariant controller that achieved global exponential position 
tracking but the vessel orientation could not be controlled. In addition, it is shown in Q that the 
underactuated ship can not be transformed into a driftless chained system; which means that the 
control techniques used for the similar problem of nonholonomic mobile robot control cannot be 
applied directly to the underactuated ship control. Accordingly, control of underactuated vessels 
in this configuration has been studied rigorously by contemporary researchers, examples of which 

are H @, 103, CD, E). 

In 0, the author showed that under discontinuous time- varying feedback, the underactuated 
vessel is strongly accessible and small-time locally controllable at any equilibrium. A discon- 
tinuous time-invariant controller was proposed which showed exponential convergence of the 
vessel towards a desired equilibrium point, under certain hypotheses imposed on the initial 
conditions. In [|TJ, a continuous periodic time- varying feedback controller was presented that 
locally exponentially stabilizes the system on the desired equilibrium point by using a global 
coordinate transformation to render the vessel's model homogenous. In a combined integrator 
backstepping and averaging approach was used for tracking control, together with the continuous 
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time-varying feedback controller for position and orientation control. This combined approach, 
later on used in [13], provides practical global exponential stability as the vessel converges to a 
neighborhoo d of the desired location or trajectory, the size of which can be chosen arbitrarily 
small. Jiang [14J used Lyapunov's direct method for global tracking under the assumption that the 
reference yaw velocity requires persistent excitation condition; therefore implying that a straight 
line trajectory could not be tracked. This drawback was overcome in |[T5l and [fT6ll . Do et al. [fT51 
proposed a Lyapunov based method and backstepping technique for stabilization and tracking 
of underactuated vessel. In this work, conditions were imposed on the trajectory to transform 
the tracking problem into dynamic positioning, circular path tracking, straight line tracking and 
parking. Borhaug et al. have proposed a control scheme for straight line following of a formation 
of marine vessels in IfTTl . 

In this paper, we address the global tracking control of underactuated vehicles, using saturated 
state feedback control OH, OH, [HQl [ED, IE21, fl23]l. Our work addresses the remaining case 
not treated in ITT3TI . i.e., the yaw angle of the tracked trajectory does not admit a limit at time goes 
to infinity. This research is therefore in the same direction as in 11241 . where the author achieved 
practical stability. Our algorithm provides asymptotic convergence to the tracked trajectory from 
any initial point. The advantage of using saturated controls is that the global asymptotic stability 
is ensured while the control inputs remain bounded, as real life actuators are all limited in 
output, see for instance [|25l . I1261 . The proposed controller has been proven to work with state 
measurements, as well as with observers in the case where all states may not be measured (cf. 
also Il25ln . 

The paper is organized as follows: the vessel model is presented in Section 2 and the control 
problem is formulated in Section 3. In Section 4, the controller is developed and the proof of 
stability is given. In Section 5, the stability of the controller is shown in presence of observation 
errors. Simulations are given in Section 6 and concluding remarks are presented in Section 7. 

II. Vessel Model 

In this section, we discuss the physical model of the marine vessel and the related assump- 
tions on physical phenomena associated with its motion. Then, a mathematical reformulation is 
presented, following variable and time-scale changes, to obtain a suitable form for control design. 

A. Physical Model 

The general 6-DOF rigid body model for surface marine vessels presented in [|2l can be reduced 
by considering surge, sway and yaw motions only, under the following assumptions [|24l . 
(Hi) Heave, roll and pitch motions induced by drift forces of wind, wave and ocean current 
are neglected. 

(H2) The inertia, added mass and hydrodynamic damping matrices are diagonal. 
The aft propeller configuration provides only the surge force x u and the yaw moment x r . The 
kinematic and dynamic equations of the vessel can therefore be written as 03J, (|24| 

{x = wcos(ty/) — vsin(i//) , y = Msin(ty/) +vcos(ty/) , 
u = -vr — au + ii, v = —cur — bv, (1) 
c 
l/A = r, r = Kuv — dr-\-T2, 

where (x,y) and i// are the coordinates and the yaw angle of the vessel in the earth- fixed frame, 
and u, v and r denote the surge, sway and yaw velocities respectively. The control inputs f\ and 
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%2 are the normalized expressions of the surge force and yaw moment, given as Ti 
1 



mi 



-T„ and 



T2 = — T r . The parameters a, b, c, d and K are positive constants that represent the mechanical 

properties of the system, namely the inertia m, > and hydrodynamic damping di, where i = 
1, 2, 3 corresponds to surge, sway and yaw motions respectively. The constants are defined as 



follows a 



d\ 



,b 



d 2 
ni2 



mi 

m,2 



, d 



d-x 



m 3 



mi — m 2 
m 3 



B. Model for control 

For control design, the system model ([[]) can be simplified by normalizing the physical 
parameters through straightforward variable and time-scale changes. For the sake of clarity, let 
us rewrite System ([T]) as follows, 




R 



-D 



u 



VI c \ V 



-rA 

r, r = Kuv — dr + %2 
where the matrices Do, Ry and A c are given as 

a 
b 



u 



+ 



(2) 



o 



R 



cos(yA) -sin(v) 
sin(yA) cos(i^) 



-1/c 




(3) 



Let us consider the following matrix D p = diag(p,cp), where p is a positive constant to be 
chosen later. Then we obtain A\ = D p l A c D p . The time scale s: = dt is introduced in System ([2]) 
as well as the linear changes of variables (u(t) / dp ,v(t) / dcp) and r(t)/d, still denoted (u(s),v(s) 
and r(s) respectively. 



After easy computations and by setting /3 : = 



K 



cp z 



r, Ti := 



pd 2 



^2 

d 2 



and D = diag(a/d,b/d), 



the dynamics of the vessel, denoted by (5), is rewritten as follows, 



(S){ 




RyDp 



-D 



r, 



u 

V 

II 

V 

f = fiuv 



rA 



ii 

v 

r+T 2 - 



(4) 



III. Problem Formulation 

The goal of this paper is tracking control of the presented underactuated marine vessel by 
controlling its position and orientation. The vessel is forced to follow a reference trajectory 
which is generated by a "virtual vessel", as follows E71 . Il28ll . 



(Sre) { 




R Vre D P 



D 



u re 

Vre 



r re Ai 



PUr e V n 



Vre 



+ 



(5) 



,rei 
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where all variables have similar meanings as in System (Q. Tracking control is achieved by using 
saturated control inputs and under the assumption that the velocities are bounded Il29ll , ||30ll . This 
assumption holds true physically as resistive drag forces increase as the velocity increases and 
therefore the latter cannot increase indefinitely if the control is bounded. These assumptions are 
also valid for the reference system and are formalized in the following manner: 

Assumption 1. There exist constraints on the control inputs and velocities such that 

1^1 1 < tl,max, 1^2 1 < T2,max, W\ < "max, |v| < Vmax, (6) 

where Ti ;inax , T2,max> "max and v max are known positive constants. The velocities u re , v re and 
the forces X\^ re and %2,re verify the same bounds as above and the reference angle \\f re does not 
converge to a finite limit as t tends towards infinity. 

The variable and time-scale change defined in the previous section requires the following 
new bounds to be defined for the new control inputs l\ and T2, denoted by Ti max and T2, m ax 

max „ T . , j. „ . .... 

respectively: X\ max = t~ an d ^2 max = — ~=r~ ■ We consider the following condition upon the 

pd z d z 

saturation limits of the control inputs, to be used later on in the control design. We use here m\ 
to denote mm{a\j2,b\). 

T 2 

o I, max „ fn\ 
CI: (5— <T2,max- ( 7 ) 

a\m\ 



max 1 H r-\ 1 

Note that this condition is always satisfied by choosing p > — - — W . Our control 

a y aimiT2,max 

objective is that (5) follows (S re ). With respect to the frame of reference of the reference trajectory 
(S re ), the error system is defined as 

e x = cos(\j/ re )(x-x re ) +sm(\j/ re )(y-y re ), e y = -sm(\j/ re )(x - x re ) +cos(v^s)(y-y«), 
e u = u-u re , e v = v-v re , e w =y-y re , e r = r-r re . 

Defining new controllers w\ and w 2 , as follows, w\ := X\ — %\ re and W2 '■= 12 — T2,re> me dynamics 
of system ([8]) becomes 



(Se){ 



+ (cos (e w ) - l)l)p ( J£ J + (R ev -Id 2 ) ( ^ ) 

















e v ) 



e r , e r = /3 (wv - u re v re ) -e r + w 2 . 



(9) 



The control objective is to force the error system (S e ) to 0, using w\ and #2. 

IV. Controller design 
We first develop the following intermediate result, concerning the bounds of u, v, r. 

Lemma 1. The variables u, v, r are bounded and satisfy 

limsup||(w,v)|| < . Tl .' max , limsup \r\ < T 2 , max + P - 1,max ■ (10) 
t^oo 2y/a\m\ t^oo 2a\m\ 
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Proof: Let us consider the positive definite function 

K,v = ^ 2 + v 2 ), 

then 

V U)V = uu + vv = u (— a\u + rv + Z\) + v (— b\v — ru) 
a\ 9 , 9 T 2 
2 2ai 
T 2 

< -m l (u 2 + v 2 ) + -^-, 

L(X\ 

X 2 

we can see that V uv < for (w 2 + v 2 ) > - — - — , which directly implies the first inequality in 

1a\m\ 

( fTOj ) (see Chapter 9 of fllQ)- Similarly, consider V,- = -r 2 then, K = rr = -r 2 - r (x 2 + /3wv) = 
— r(r + T2 + /3wv), V,- < for \r\ > | T2 1 H— |/3wv|, from which we derive the second inequality in 



era). ■ 

Remark 1. As the reference trajectory system is similar to the vessel model, it can be shown 
that the limits defined in Lemma [7J are valid for well. 

We define a new control variable W2'-=P (uv — u re v re f) + W2. As the upper bounds of u, v, u re 
and v re f are known according to Lemma [T] and Remark [Tj we obtain 

T 2 

limsup/3 \uv — u re v re f\ < /3 — — . (11) 
r^oo 1 1 a\m\ 



X 2 max 

If |w 2 | < U2, for a positive constant U2, one must have U2 + /3 — ! — < T2 m(U , which is guaranteed 

a\m\ 

by Condition CI. With these preliminaries established, we will now proceed to fulfill the control 
objective by using the bounded controls w\ and W2- Let a(.) be a standard saturation function, 

i.e., o(t) = — — . The main result of the paper is given next. 

max(l, \t\) 

Theorem 1. If Assumption [7J and Condition CI are fulfilled, then for an appropriate choice of 
constants U\, p, A/, U2, k\, &2, M satisfying: 



a x > Ui +p, Ui > 



P, c/ 2 >o, 



min^ai, ^ 

M > 0, k\ > k 2 - 1 > 0, ai +«§ = /ip, &i = jUcp, 
then the following controller ensures global asymptotic stability of the tracking error system (S e ): 



^1 ^2 — 1 

t>2 U2 



(12) 



Proof of Theorem HJ 
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We first consider the errors ey and e r and take k\ , ki > large in the control input W2 defined 

previously. The dynamics of e v and e r in (S e ) ey = e r and e r = —e r — U2O ( jj-e v + ^ 2 - e r ) . 

\U2 U2 J 

Lemma 2. t/2 > and k\ > Ic2 — 1 > 0, then after a sufficiently large time, the saturated control 
operates in its linear region and the errors ey and e r converge to zero exponentially. 

Proof: Let us consider the Lyapunov function V, 



a 9 
V = -e 2 r +S 



k\ k2 — l 



(13) 



where S(cj) := / o(s)ds, and a 
Jo 



ki — k2 + 1 k\ k2 — 1 

~ > 0. Set z := — ey-\ e r . Then, one has 

Ui U2 U2 



V = ae r e r + o(z)z = -cce 2 - (k 2 - l)o 2 {z). 



Then V < for (ey,e r ) 7^ (0,0); and after a finite time we obtain 



h 
U 2 



k 2 -\ 



(14) 



< 1. The 



dynamics of ey and e r becomes linear, i.e., ey = e r and e r = —k\ey — k2e r . As £i,&2 > 0, (ey,e r ) 
converges exponentially to zero. ■ 
Lemma [2] shows that the errors ey and e r converge to zero under the control W2- We will 
now consider the errors e u and e v . We choose the constants /i and t. such that a\+£> = jlp and 
b\ = jlcp, implying that t, =b\/c — a\ and \i > 0. 



Lemma 3. Consider the dynamics of e u and e v given in Equation (|9]). If U \ and p are chosen as 

b\ 

( *r ^ 

-p, then the control w\ := —U\0 



a\> U\+ p and U\ > 



a\ — 



mm ( a\ , -f 



-p<7i(.), with CTi(.) 



to Z?e chosen later, ensures that e u and e v satisfy the following inequalities: 



limsup||(e M ,e v )|| < -^=== : 



wzY/z a = inf I a\ + E 
t>0 1 




> 0, m2 := min(a/2,fci). (15) 



Proof: Notice that a > since it is trivially the case if % > 0, and otherwise, <5 > «i + 
£ = — . Then, one has e„e M + e v e v = -a^ — b\e% + e r {e u Vre - e v u re ) + e u W\. By applying the 

control wi, we get e B e B + e v e v < -a\e\ - b\e\, - U\e u o I J - pe u O\(.) + C |e r | ^e 2 + e 2 , 

where Co := iimax + Vmax- According to Lemma [2j e r tends to zero, i.e., for large t, one has 
e u e u + e v v < — [a\ + %s(t)] e 2 — b\e 2 — pe u O\ (.). Then ([15]) results from the following inequality 



e u e u + e v v < -m 2 {e 2 u + e 2 ,) + 



(16) 



a 



Lemma 



one gets that lim sup 



proves the convergence of e u and e v to a neighborhood of zero. Since 



~> • ( b ^\ 

a > mm I a\, — I . 



< 1, and the controller exits saturation in finite time and enter its linear 
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region of operation. We get o 



Ui 



, and the dynamics of e u and e v become 



e u = -npe u + r re e v - p<J\{.) + e r v re - e r e v , e v = -[icpe v -r re e u - e r u re + e r e u . 
Define W = (W U W 2 ) T := (e x ,e y ) T + (e u ,e v ) T /^l. Then one has the following result. 



Lemma 4. Let W defined previously and the controller w\ given in ( fT2[ ) with Oi(.) = o(MW\), 
where M is an arbitrary positive constant. Then W tends to a finite limit W = (0, W2Y ■ 



Proof: The dynamics of W can be expressed as 



W = -r re A i W + —A l 



e v 



— r. 



+ sin (e w )AiD p 
P 



o(MWi 



U re 
Vre 
1 





+ 



+ 0(\e r \(l + \\(e u ,e, 



In order to find the limsup of ||W||, we calculate 



W T W 



sin (e ¥ ) W T A l D p 



+ W 1 



01 e 



w I vl 



e v 



+ 0(|^|) 



Wio{MWi), 



0(\\W\\.\\e v ,e r \\)--W l o(MWi] 



which implies that \W T W\ + — WiOi(MWi) ^ ||W||0 (||e v ,e r ||). One deduces first that the time 

derivative of ||W|| is integrable over E + and thus W admits a limit as t tends to infinity. Therefore, 
the right-hand side of the previous inequality is integrable over E + implying the same conclusion 
for W\o(MW\). As both W\ and W\ are bounded, then according to Barbalat's Lemma, Wi — > 
as t —7- 00. Consequently W2 tends towards a finite value W2 as t tends to infinity. ■ 
Lemma [4] permits us to further improve the result of Lemma [3j as follows. 

Lemma 5. If Lemmas [5] and |4] hold true, then e u and e v converge to zero asymptotically. 

Proof: From Lemma [3] and setting G(e u ,e v ) := (e„ + el)/2, Equation ( p~6| ) is rewritten as 



G + 2m2G < 



One concludes using Barbalat's Lemma. 



So far, we have established that the errors e^, e r , e u and e v converge to zero. From Lemmas 
|4] and [5j we deduce that if W\ — > and e u — > 0, then e x will converge asymptotically to zero as 
well. We next address the convergence of the remaining error variable, ey. 



Lemma 6. If Assumption \ljis satisfied, then W2 = and e y converges asymptotically to zero. 

Proof: From Equation ( fT7| ) in Lemma |4} the dynamics of W can be expressed as follows, 
W = -r re AiW + O (\e w \ ,\e r \ ,W\0 (MW{)). We define the new variable W as follows, W : = 
R\lf re W, and the dynamics of W is given by 

W = W re AiR Wre W - r^R^W +R ¥re O(\e v \ , \e r \ , W x a (MW t )) =0(\e ¥ \, \e r \ , Wio{MW x )) . 

Since and e r converge exponentially to zero and W\o(MW\) is integrable over IR + , then |[w|| 
is also integrable over E + , which means that W converges to a finite limit W l . Then, one gets that 
— VV2sin(yA re ) and Videos (yf re ) tend to W[ and W{ respectively, as t tends to infinity. If W2 ^ 0, 
then one easily shows that \\f re converges to a finite limit by considering whether (W re ) 2 = or 
not. That contradicts Assumption 2 and W must converge asymptotically to as well as e y . ■ 
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It should be noted that the controller presented in Theorem 1 has been designed under the 
assumption that all state variables are known. In the next section, the study is extended to the 
case where only the position and orientation states of the vessel are available and the velocities 
need to be observed. 



V. Tracking without velocity measurement 

In practical cases, only position and orientation feedbacks are available for navigation. Therefore 
the only available states of the vessel are x, y, ty/ along with the the complete coordinate state 
set of the virtual vessel to be followed. For such output feedback systems, the variables u,v,r 
need to be observed. In this section, we will show that the controller presented in Theorem [T] 
is applicable in this case and the use of observation instead of measurement does not affect the 
stability. We suppose that the velocities are obtained through an observer such as that presented 
by Fossen and Strand in ll32ll . or a robust differentiator such as that presented in [33]. In both 
cases, observation errors converge exponentially to zero. It can be noted that, when we use 
a differentiator, the estimated values (w,v,r) of (w,v,r), can be determined according to the 
following equation (u,v) T = Dp R-y,(x,y) T and r=\j/, where (x,y,\jf) are the estimated values 
of (x,y,yf) respectively. 

Let us follow the same steps used in the demonstration of stability of system (S e ) with velocity 
measurement. The observation error related to the velocity are define as below: f u = u — u, 
f v = v — v and f r = r — f. As the references are common, the observation errors can be described in 
terms of trajectory pursuit errors as /„ = e u — e u , f v = e v — e v , and f r = e r — e r , where, e u = u — u re , 
e v = v — v re and e r = r — r re . We note that the variable x,y, y/ are measured and the related 
observation errors are null. The problem is transformed to demonstrate the stability of the error 
system (S e ) under control laws w\ and #2, which are now based on the observed values. 
As in the previous case, we define W2 '■= /3 (uv — u re v re ) + W2- Then the result of this section can 
be stated as the following theorem. 

Theorem 2. If Assumption [7J and Condition CI are fulfilled, then for an appropriate choice 
of constants U\, p, M, U2, k\, Ic2, jU, the following controller ensures global asymptotic 
stability of the tracking error system (S e ): 

m = -Ukj [^j -pa \M(e x + ^e u yj ,w 2 = -U 2 o ^^^) > ( 17 ) 

if in addition the observer errors f u , f v and f r converge asymptotically to zero and are integrable 
over IR + (i.e., the integrals of their norms over R+ are finite). 

Remark 2. The choice of constants U\, p, M, U2, k\, &2, M remain the same as in the case 
of Theorem [7] therefore their expressions and conditions will not be repeated in this section. 

Proof: The proof of Theorem [2] is largely based upon the proof of Theorem [TJ and is 
developed similarly. We first consider the dynamics of error variables e^ and e r : 

e ¥ = e ri e r = -e r -U 2 a(^e v + 1 e r + f\ (?) j + f 2 (t) , (18) 

k 2 - 1 

where f\(t) = — - f r and f 2 (t) =fi(uv-uv). 

Ll2 

Lemma 7. If f\{t) and f 2 (t) converge to zero asymptotically, then for some large positive 
constants k\ and k2, System (fT8l) is globally asymptotically stable. 
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. . k\ k2 — 1 

Proof: Consider the Lyapunov function V defined in ([13]). If z := jr e V^ — Ti — one § ets 

U 2 L>2 



, k\ k2 — 1 

V = ae r e r + o | — e v + — — e, 
U 2 U 2 



U 2 



k 2 -\ 
U 2 



ae r {-e r - U 2 o{z + fi) + f 2 ) + o(z) 



k\ ko — 1 / 

^er+^j- {-e r -U 2 o{z + fi)+h) 



(19) 



. . a 2 + b 2 k\ - k 2 + 1 
Using the inequality, \ab\ < — - — , and taking (1U2 = > 0, we get 



U 2 

V < -^e 2 r -(k 2 -l)a(z)a(z + fi) + aU 2 e r (a(z)-a(z + f l )) + ^fi + k2 ' 1 



a 

' — 1 

2 



U 2 



\fl\- 



Using the inequalities \a(z) — o(z+fi)\ < \fi\ and o(z)o(z + fi) > o 2 ' (z) — |/i|, one has 



v<- rr 



a 2 (k2-l)a 2 (z) + (k2-l)\fi\ + aU 2 \e r \\f l \ + ^f 2 + k2 1 



U 2 



1/2 



After a sufficiently large time interval T, it is assured that \U 2 fi \ < - V? > T, and 

6 



a 



V < --e 2 - (k 2 - l)o 2 (z) + O (ff, \f 2 \,fi) . 
From here, we obtain that limsup \e r \ = limsup |<r(z)| = 0. 



(20) 



Following the same steps as used in the previous section, we now demonstrate the convergence 
of the error variables (e u ,e v ,e x ,e y ) of System (S e ). 



Lemma 8. Consider the dynamics of e u and e v presented in Equation 



-Uio 



•poi(0 



Then, the control 

(21) 



ensures that e u and e v are bounded and again verify the estimate (15). 



Proof: The argument exactly follows the line of the argument of Lemma [3] with the controller 
-U\0 f -77^+/ J — POi(.), where / := — — /„ converges to zero asymp- 
, and by using the inequality xo(x + f) > xo(x) — \xo(f)\. 



w\ written as w\ 



totically and is integrable over 



Similarly, the proof of convergence of W\ to zero and W 2 to a finite limit, presented in Lemma 
[6| holds true and one concludes by showing that the limit of W 2 is zero as well. ■ 



VI. Simulations 

The performance of the presented controller is illustrated by simulation. We apply the controller 
on an example of a monohull vessel, as considered in [15]. The length of this vessel is 32 m, 
and a mass of 118 x 10 3 kg. The parameters of the damping matrices as given as follow: 



d\ 
m\ 



215 xl0 2 £gs-\ d 2 = 91xl0 3 Kgs- 1 , d 3 = S02xl0 4 Kgm 2 s 



2„-l 



120x10-^, m 2 



172.9 x 10 j £g, m 3 



636 x l0 5 Kgm 2 . 



Based on these physical parameters, we find the parameters of System ([T]) as 
a = 0.179, b = 0.561, c = 0.694, j8 = 0.126, K = 8.32 x KT 4 . 



(22) 



(23) 
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Then, the parameters of the controller and the normalized system (5) are given by 

a x = 1.421, bi= 4.449, d = 0.126, £i = 10, fc 2 = 10, t/ 2 = 0.1, t/i = y , P = y , M = 0.1. 

The reference trajectory is generated by considering the surge force and the yaw moment as 
constants Ti >re = 10 and %2,re = 0.05 with the initial values 

(x re (0),y re (0), V^e(0),w,- e (0),v re (0),r re (0)) = (0 m,0 m,0 rad, 0ms 1 ,0 ms~ l ,0 rads~ 1 )). 
The initial conditions of the vessel are as follows: 

(x(0),y re (0), y(0), w(0), v(0), r(0)) = ^50 m,-150 m, | rad, 50ms~\0 ms~ l ,0 rads^ . 

The reference trajectory and the vessel are shown in a 2D coordinate plane in Figure [TjThe vessel 
converges to the reference trajectory asymptotically and similarly for the position errors graph 
in Figure [2j The orientation error and its derivative also converge to zero, as seen in Figure |4j 
The convergence of e u and e v is shown in Figure [3} Figures [5] and [6] show the control signals Ti 
and T2 respectively and the controllers are clearly bounded. This is an essential property in real 
systems, where the control signals are constrained. 

VII. Conclusion 

In this paper, we have addressed the problem of tracking of an underactuated surface vessel 
with only surge force an yaw moment. The proposed controller ensures global asymptotic tracking 
of the vessel, following a reference trajectory modeled by a virtual vessel. It is also shown that 
the stability of the system is not affected if the state measurements are replaced by observers. The 
using of saturated inputs is essential as in real life the actuators have limited output. Simulation 
results illustrate the performance of the controller. 




Fig. 1. Reference trajectory and the vessel 



Fig. 2. Errors e x and e 
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time time 

Fig. 5. Control Ti Fig. 6. Control T2 
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